Predicting scattering properties of ultracold atoms: adiabatic accumulated phase 

method and mass scaling 
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Ultracold atoms are increasingly used for high precision experiments that can be utilized to 
extract accurate scattering properties. This calls for a stronger need to improve on the accuracy 
of interatomic potentials, and in particular the usually rather inaccurate inner-range potentials. A 
boundary condition for this inner range can be conveniently given via the accumulated phase method. 
However, in this approach one should satisfy two conditions, which are in principle conflicting, and 
the validity of these approximations comes under stress when higher precision is required. We show 
that a better compromise between the two is possible by allowing for an adiabatic change of the 
hyperfine mixing of singlet and triplet states for interatomic distances smaller than the separation 
radius. A mass scaling approach to relate accumulated phase parameters in a combined analysis 
of isotopically related atom pairs is described in detail and its accuracy is estimated, taking into 
account both Born-Oppenheimer and WKB breakdown. We demonstrate how numbers of singlet 
and triplet bound states follow from the mass scaling. 

PACS numbers: 34.20.Cf,34.50.-s,67.85.-d,67.85.Fg 



I. INTRODUCTION 

In 1976 Stwallcy H suggested the existence of mag- 
netically induced Feshbach resonances in the scattering 
of cold hydrogen atoms. He pointed out that the spe- 
cific magnetic field strengths where they occur should be 
avoided to achieve a stable cryogcnically cooled H gas, 
in view of an enhanced decay at resonance. In 1992 one 
of the present authors (BJV) and co-workers pointed to 
a positive aspect of such Feshbach resonances they 
allow for an easy control of the interaction strength be- 
tween ultracold atoms, i.e., atoms in the energy range 
where their interaction is limited to s-waves. In such 
circumstances, the interaction strength is characterized 
by the s-wave scattering length a. With a Feshbach reso- 
nance, the interactions can be tuned from weak to strong 
and from attractive to repulsive by simply changing an 
externally applied magnetic field. 

Since then these resonances have become an indispens- 
able tool in many successful attempts to control the in- 
teratomic interaction, to form ultracold molecules by as- 
sociating atoms, and to create a superfluid Fermi gas. 
Feshbach resonances allow experiments with ultracold 
atoms access to a multitude of the most diverse many- 
body phenomena 0] . Systematic theoretical work to de- 
termine resonant field strengths and scattering lengths 
for almost all stable alkali atoms started immediately 
after 1992 0, i, H, 0] and played a crucial role in the 
first realizations of Bose-Einstein condensation BEG in 
1995 [1, i, [H (we will use the term alkali atom in ref- 
erence to an alkali metal). An example is presented in 
Sec.|ll]in connection with the first determinations of scat- 
tering lengths. In recent years many experiments have 
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opened the field of ultracold gases with mixed atomic 
species, where Feshbach resonances continue to be an in- 
dispensable tool. 

A description of cold collisions between ground-state 
atoms (and also weakly bound states) requires highly 
accurate central interaction potentials. Except for the 
lightest elements (H and Li), ah initio potentials do not 
possess the required accuracy at short range. The slight- 
est change of a potential in that range can easily turn 
a positive into a negative scattering length, information 
which is crucial for instance to predict the stability of a 
BEG. 

A way to account for that is to summarize the "his- 
tory" of the collision for interatomic distances r smaller 
than a separation radius rg by means of a boundary con- 
dition on the wave function at 7'o , and to determine that 
condition from a restricted set of available experimen- 
tal data [1, 0, @1 ■ The basic philosophy of this approach 
is to give up the goal of extracting the detailed short- 
range potential as a whole from experiment in favor of 
a boundary condition with only a few parameters. The 
boundary condition takes the form of a radial phase of 
the zero-energy wave function accumulated in the inter- 
val r < To in either the singlet or the triplet channel, and 
its energy and angular momentum derivatives. This pre- 
supposes pure singlet and triplet wave functions, which 
is justified for small interatomic distances where the sin- 
glet and triplet states are far enough apart in energy to 
neglect hyperfine mixing. 

Over the years the accuracy of the description of scat- 
tering properties obtained with this method has shown a 
dramatic improvement, keeping pace with the accuracy 
of the measurements. In this paper we describe an exten- 
sion of the accumulated phase method, the adiabatic ac- 
cumulated phase method, presented briefly in a previous 
publication [ll| , in which the hyperfine interaction is not 
completely neglected, but taken into account adiabati- 
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cally for r < tq. It enables one to describe the interac- 
tion and scattering of (ultra)cold atoms to unprecedented 
precision. It is also unparalleled in comprehensiveness: 
it allows the prediction of a large and varied set of ex- 
perimental data, once the accumulated phase parameters 
have been determined from a restricted set of experimen- 
tal data. We thus refrain from extracting the short range, 
deep part of the central interaction potentials, singlet and 
triplet, in favor of predicting new scattering properties. 
As such, it sets an example for future " state of the art" 
applications to other atoms, in particular to interactions 
between unlike atomic species. It should also be noted 
that predictive power of the accumulated phase method 
extends to both scattering and bound states, because the 
corresponding wave functions do not differ essentially at 



As a background for the later exposition we start in 
Sec. lUwith a summary of the accumulated phase method 
as we introduced and applied it previously. This sec- 
tion also serves to introduce our notation for the various 
terms in the Hamiltonian. As a last step in that section 
we describe our approach to include the spin-spin inter- 
action, in particular its second-order spin-orbit part, in 
the accumulated phase method. The adiabatic extension 
of the latter method is presented in Sec. lIIII We compare 
the effectiveness of the new adiabatic accumulated phase 
method to that of the conventional method. Sec. lIVI deals 
with the subject of mass scaling. It is based on the Born- 
Oppenheimer and WKB approximations. In this connec- 
tion it should be emphasized that the concept of mass 
scaling as studied here is basically different from that in 
other studies of cold atom scattering (see, e.g. Ref. [T3|). 
in that we apply it to a restricted range of interatomic 
distances thus avoiding the range r > rg, in part of 
which the central potentials become too shallow to allow 
for an accurate scaling close to the dissociation energy. 
In Sec. |V] we discuss the accuracy of the mass scaling, 
taking into account both Born-Oppenheimer and WKB 
breakdown, thus showing mass scaling to be a promising 
method to relate accumulated phase parameters for dif- 
ferent isotopes of the same element. Sec. IVII compares 
the conventional and adiabatic accumulated phase meth- 
ods and discusses the dependence of our predictions. 
In Sec. IVIII we demonstrate our method of determining 
the numbers of pure singlet and triplet bound states of 
the various atomic species involved in the mass scaling. 
To our knowledge this is the only available method to 
extract numbers of bound states without knowing the 
potentials at short range. We illustrate it by means of 
the example of the ®^Rb and ^^Rb atoms. A summary 
and outlook are presented in Sec. IVIIII 



II. INTERACTIONS AND ACCUMULATED 
PHASE METHOD 

Single particle Hamiltonian 

We consider the electronic ground state of an alkali 
atom. The valence electron has spin s=i while the nu- 
cleus has spin z, which in particular equals | for ^^Rb and 
I for ®^Rb. Note that lower case characters are used to 
indicate single atom properties while we reserve capitals 
for two-atom systems. In total there are thus 2(2i -I- 1) 
possible 'ground states' for an alkali atom. This degen- 
eracy is lifted by interactions both within the atom and 
with external fields. To begin with, the nuclear spin in- 
teracts with the valence electron spin via the hyperfine 
interaction for an atom j: 

y^'-^^r^^^ (1) 

with a^^ the hyperfine constant (dimension of energy), 
leading to the hyperfine splitting according to / = z ± ^ 

(/ = s + i) with each /-state (2/-f l)-fold degenerate. In 
turn, the remaining degeneracy is lifted when atom j is 
placed in a magnetic field B, giving rise to a spin Zeeman 
term 

Vj^ ={le^l^l-lnJj)-B, (2) 

where 7ej- and 7„j are the electronic and nuclear gyro- 
magnetic ratios. The behavior of the valence electron in 
alkali atoms is influenced by the electrons filling the inner 
shells causing the gyromagnetic ratio ^e,j to be slightly 
different from that of a free electron. Eqs. H]) and ^ 
lead to the familiar graphs for the energy of the atomic 
hyperfine states as a function of the magnetic field Bz 
(see Fig. [1]). The single atom hyperfine states are com- 
monly denoted by |/, my) even though / is only a good 
quantum number for _B = 0. 

Two particle Hamiltonian 

We consider two like alkali atoms. They experience in 
the first place a mutual central interaction, that can be 
written as 

y-"(r) = Vs{r)Ps + VT{r)PT, (3) 

with Ps,T projection operators on the spin singlet and 
triplet subspaces and r the interatomic separation. The 
singlet and triplet potentials differ by twice the exchange 
energy V'^'^'^^{r) and are at large distances given by 

Vs^T = V^'^'P - (-1)^ V'"'"'''. (4) 
The dispersion energy y^'^P(r) is described by 
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FIG. 1: Hyperfine diagram of ground-state *^Rb. The energy 
curves of the hyperfine states are labeled by their B=0 quan- 
tum numbers. A similar diagram with i = 5/2 is valid for the 
**^Rb atom. 



with the dispersion coefficients C„. An analytic expres- 
sion for the exchange energy in Eq. (j4|) has been derived 
by Smirnov and Chibisov [l3| for r values where the over- 
lap of the electron clouds is sufficiently small: 



ycxch ^ i 



(6) 



In this equation J and k are positive constants with k? /2 
the atomic ionization energy; r, J, and k = 0.554.... in 
atomic units. The most recent value for J was given by 
Hadinger et al. , who made use of Ref. . 

Leaving out the center of mass kinetic energy and in- 
cluding the above interactions the total effective Hamilto- 
nian for two colliding ground-state alkali atoms becomes 



H 



2n 



(7) 



in which the first term represents the kinetic energy with 
fi the reduced mass and p the relative momentum oper- 
ator. 

The hyperfine term can be written as the sum of 
two parts with different symmetry with respect to in- 
terchange of the electronic or nuclear spins, 
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The convenience of this splitting arises from the fact that 
V^^~^ is diagonal in S, whereas V^^~ , being antisymmet- 



ric in si and 5*2, is the part coupling singlet and triplet 
states. 

For the interactions mentioned up to now the total 
Hamiltonian H is invariant under independent rotations 
of the spin system and the orbital system around the 
axis through the overall center of mass parallel to the 
magnetic field. Therefore the projection of the total spin 
angular momentum fi+f2 = F and the orbital angu- 
lar momentum I along this axis are separately conserved 
during the collision. Since V'^™ depends on r only and 
not on f = f/r, I is even conserved as a 3D vector. As a 
consequence, mp and the rotational quantum numbers I 
and mi are good quantum numbers. 

Two other interactions influence the two-atom system, 
which are much weaker than the above-mentioned effects, 
but nevertheless can play a significant role in interpreting 
cold atom experiments. The first one is a direct interac- 
tion between the spins of the valence electrons via their 
magnetic dipolc moment. It is given by 



Ml ■ - 3(/ii ■ f)(^2 ■ r) 



(9) 



with Mo = 47r • lO^^Hm^^ and jlj the electronic magnetic 
dipole moment of atom j. We leave out the much weaker 
magnetic dipole interactions in which the nuclear mag- 
netic moments are involved. Another interaction with 
similar structure arises as a second order effect of the 
spin-orbit interactions of the valence electrons [3 via 
an intermediate coupling to electronically excited states. 
The well-known electronic (first-order) spin-orbit inter- 
actions do not play a role for the orbital s-states of the 
valence electrons. 

In total, we thus have a spin-spin interaction V^^ be- 
tween the valence electrons, consisting of two parts: 



(10) 



a magnetic dipole part and a part arising from the spin- 
orbit interaction in second order. The dipole part, when 
expressed in the spin vectors si is given by 



M07e 



[Si ■ 32 - 3(si • f)(s2 • f)] 



(11) 



The part (1/^**)^° has effectively the same spin-angle 
structure (the factor between square brackets), but is 
multiplied by a different radial factor. This factor has 
been calculated via an ah initio electronic structure cal- 
culation by Mies et al. and can be approximated as 
an exponentially decaying effective form /(r) for increas- 
ing r: 



1 



fir) = --^C-'"^^ 



(12) 



with a the fine structure constant and C"", (3 and r^o fit 
parameters to the ab initio results. In the following we 
will show that its effect on two-atom bound or scattering 
states can effectively be accounted for via one parameter 
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only [ii,[Til, which has the form of an integral of f{r). 
The total V^^ apparently has the spin-angle structure of 
a scalar product of two irreducible spherical tensors of 
rank 2 ^ (Sec. 31): 

2 

{si,S2)^- ir,r)2 = ^ (-1)" (5*1,32)2^ (f,f)2_^ 

2 

^ ^(-l)-(.-i,.-2)2.y2-.(r)(13) 



As a consequence, it is invariant under the simultaneous 
3D rotations of the orbital and spin degrees of freedom, 
thus conserving the total two-atom angular momentum. 
On the other hand, it is not invariant under independent 
orbital and spin rotations. It therefore obeys triangle 
type S and I selection rules for a second rank tensor: it 
couples only spin triplet states and it couples for instance 
the / = and 2 relative partial waves of the two atom 
system. 

As mentioned before, both parts of the spin-spin in- 
teraction are relatively weak. In most applications they 
can safely be neglected. That these interactions cannot 
always be neglected is illustrated by the observation of 
Feshbach resonances in ^'^^Cs [2l| and ®^Rb [2^ exper- 
iments in which colliding ultracold atoms, approaching 
each other in an s-wave, resonate with an 1=2 or even Z=4 
(quasi-)bound state coupled via the spin-spin interaction 
to the entrance channel (in the l^i case this interaction 
is needed twice: s-wave d-wave f/-wave). 



Accumulated phase method 

As a background for the later exposition of its new adi- 
abatic variant and the mass scaling procedure, it is useful 
first to recapitulate the basic features of the conventional 
accumulated phase method. The separation radius tq 
where we impose the boundary condition is chosen so as 
to fulfill three conditions: 

1. rg should be so small that in the range r < vq the 
lowest S = and S = 1 two-atom electron states 
(see Fig. [2] and Fig. [3] for a pair of Rb atoms) are 
sufhciently far apart in energy for the singlet-triplet 
coupling due to V^^^~ to be negligible. This makes 
it possible to formulate the boundary condition in 
terms of pure singlet and triplet waves. 

2. On the other hand tq has to be so large that the 
singlet and triplet potentials for atomic distances 
r > rg can be accurately described by their asymp- 
totic form V'^'^'P T V'''"'^^ according to Eqs. © and 
(|6]), with a small number of imknown parameters. 

3. The value of tq, as well as both the energy E rela- 
tive to threshold and the angular momentum / val- 
ues playing a significant role in the experimental 



data should be small enough that a rapidly con- 
verging expansion of the S = and S = 1 phases 
in powers of E' and l{l + l) is possible, thus also con- 
taining a small number of unknown parameters. 

In view of the possibility that these conditions are con- 
tradictory, it is far from obvious that a suitable tq value 
can be found. In the first half of the nineties when three 
U.S. experimental groups attempted to create a BEG in 
an alkali atomic gas, it was possible to predict the signs 
and (in some cases rough) magnitudes of the scattering 
lengths for almost all alkali species, determining the sta- 
bihty (a > 0) or instability (a < 0) of a large BEG. This 
essential information could be obtained with the accu- 
mulated phase method using a value 19 and even 20ao 
for ro (ao = Bohr radius = 0.5291772 xlO"i°m). For 
example, a predicted negative a for *^Rb and a positive 
a for *^Rb atoms (both spin-stretched) led Wieman, 
Gornell and co-workers in 1995 to switch from *^Rb to 
^''Rb in their experiment, leading to the first successful 
realization of BEG in an ultracold atomic gas P . 

The concept of an accumulated phase was originally 
introduced in the spirit of the WKB approximation as the 
local phase of a rapidly oscillating radial wave function at 
rg. Its value (j)s{E, I) and (priE, I) for each of the singlet 
and triplet wave functions is defined by 



V'(ro) = A 



sm[cl){E,l)] 



(14) 



and its radial derivative, with up to a constant the singlet 
or triplet accumulated phase 



(t){E,l) = / k{r)dr. 



(15) 



Here k (r) is the local radial wave number for the channel 
involved: 



(r) 
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with /i the reduced mass and V{r) the singlet or triplet 
potential. With respect to condition (3) earlier in this 
section we repeat that for (ultra)cold colliding atoms 
(T < l/xK) and near-dissociation bound states we are 
most often considering, E is close to (compared to the 
depth of the potential at tq) and I is at most 4. As shown 
in Fig. [3]for Rb atoms, the small E and I ranges then al- 
low a first order Taylor expansion for (t>{E, I) according 
to: 



^(EJ) = 0(0,0) 



90 
'dE^ 



E 



d[l{l + l)] 



1(1 + I) 



4:" + Ecb^ + l{l + 



(17) 



The generally fractional s-wave vibrational quantum 
numbers at dissociation, vjjs and vut-, are essentially 
equivalent to the zero-order Taylor terms. They provide 
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FIG. 2: (Color online) Singlet {S = 0) and triplet (S=l) po- 
tentials for a pair of rubidium atoms in the electronic ground 
state. 

for more direct physical insight, however, being a mea- 
sure of how close the last bound or the first unbound 
two-atom state is to the dissociation threshold. Their 
fractional values arc defined via interpolation between 
successive infinite values of the scattering length making 
use of the radial phase in the deepest part of the potential 
i: 

.^(mod = (18) 

TT 

where (j>'^{a = oo) would be consistent with an infinite 
value of the scattering length, i.e. a potential which has 
a bound state at the dissociation threshold. The energy- 
derivatives correspond to the classical sojourn time 

T^o/ = 2hd(j)/dE (19) 

of the atoms in the distance range r < tq for I = and 
energies close to threshold. The l{l + 1) derivatives are a 
measure for the influence of the centrifugal force in the 
rotating two-atom system. 

It is very convenient and intuitively appealing to define 
the boundary condition in the above way. The validity 
of the WKB approximation, however, is not strictly nec- 
essary, since the phase (j){E, I) can be defined in terms 
of a logarithmic derivative. For r > tq there is a cou- 
pling region where the exchange interaction is of similar 
magnitude as the hyperfine and Zeeman energies, as in- 
dicated in Figs. [2] and [3] for the Rb atoms. For larger 
interatomic distances where V'^^'^^ has further decreased 
the two-particle hyperfine states form a good basis. 

An advantage of the accumulated phase method com- 
pared to alternatives [2^, iM, 23] is that the above set of 
phase parameters can be systematically extended by tak- 
ing more terms in the expansion (|17p into account. We 
also point to the difference with Multichannel Quantum 
Defect Theory (MQDT) methods in general: in our case 
the scattering channels are still coupled by the exchange 




FIG. 3: The S — ^ S — 1 energy splitting of two ground- 
state rubidium atoms (equal to 21/°'"^'^) versus interatomic 
separation. The hyperfine energies for the isotopes *^Rb and 
*^Rb are indicated for comparison. 

interaction in part of the exterior region r > ro where 
ycxch jg similar magnitude as the hyperfine energy, as 
indicated in Fig. [5] and [31 

Inclusion of second-order spin-orbit interaction 

In the above discussion we have neglected the spin- 
spin interactions (!/'''')'' and (y**'^)''°, which are respon- 
sible for corrections to the scattering behavior and lead 
to decay. The magnetic dipole interaction (V^^)^ con- 
tributes mainly due to its long range behavior. In 
the range r < tq the indirect spin-spin interaction (1/^**)^° 
can be very strong, exceeding the dipole interaction by 
several orders of magnitude, but it is rather weak beyond 
tq. The influence of (l/****)*^" can be taken into account in 
a model-independent way via one additional parameter. 
This can be seen as follows. In the first place we note that 
we can account for the (y'')''° mixing of the triplet chan- 
nels by means of a local S'-matrix [2^ , S_{r) , that will be 
part of the boundary condition and specifies the ratio of 
the outgoing and incoming parts of the total wave func- 
tion at r. In the vicinity of r = rp, classically accessible 
so that the local channel wave numbers are real and pos- 
itive, the radial solutions without spin-spin interaction 
are given by 

with 

4>i ^ (t)s/T {Etot - £i,h) , (21) 

where the channels i differ from each other by their sin- 
glet/triplet character, their I values or their internal ener- 
gies Ei. To formulate a local 5*- matrix at tq we introduce 
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FIG. 4: (Color online) Part A illustrates the behavior of 
the wave function's phase near ro = 16ao for three differ- 
ent energies. A comparison between the true accumulated 
phase (dots) and a first order approximation (solid lines for 
triplet, dashed lines for singlet) is shown in part B. As a 
function of E and l{l + 1) the graph shows the difference 
in accumulated phase (f>{E, Z) at r = ro as compared to the 
E = 0, I ^ situation: A(^(_E,0) = (j){E,0) - (l>{0,0) and 
A(j}{0,l) — <j){0,l) — 0(0,0), respectively. The horizontal ar- 
row indicates the typical E and / ranges for which we apply 
the first order approximation. Typical rubidium potentials 
are used for this calculation. Note that for clarity the energy 
intervals for the wave functions in part A exceed the energies 
occurring in practice by far. 



a complementary solution 



cos 



G.(r)=- 



(22) 



satisfying the Wronskian condition VF[Fi,Gi] = FiG[ — 
Fid = 1. We calculate the influence of (1/*^^)*^° by means 
of a solution of the coupled Schrodinger equation for r < 
ro using Eq. (|12p for the form factor or any alternative. 
Near ro we expand that solution in F and G functions 
according to: 



u{r) = F{r) + G{r)C (r « ro), 



(23) 





10 


- — . 




CO 




CO 


10 


E 






10 








10 












10 


O 






10 























^^^^ ^"^--0(51^155=2 11 






^""--C/CMtes=1 U 















-0.4 



-0.2 





Vdt 



0.2 



0.4 



with the non-vanishing elements of the diagonal F- and 
G-matriccs corresponding to the F and G functions. This 



FIG. 5: Rate coefficient G for the decay of Cs atoms in the 
|4, 4) hyperfine state, as a function of the potential param- 
eter vnT for several strengths of (V")^". The experimental 
boundaries for G are indicated by two horizontal dashed lines, 
the theoretical range for vdt [ig . [l9l | by two vertical dashed 
lines. 



defines a local generally non-diagonal C-matrix that ac- 
counts for the (y^*')'^° interaction. It is independent of the 
accumulated phases and depends only on the envelope of 
the rapidly oscillating F functions. Note that we could 
have used complex ingoing and outgoing exponentials as 
basis functions instead of cosine and sine functions. The 
resulting complex S_{r = ro)-matrix has a simple relation 
with C_. 

It turns out that this C-matrix can to very good ap- 
proximation also be obtained without solving a coupled 
channels equation in the range r < ro: a constant times 
the matrix representation of the spin-angular part p^ . 
independent of the radial form factor /(r), gives the 
same result. This can be understood as follows. As 
calculations show, C depends on the spin-spin interac- 
tion only in first order. A perturbation argument leads 
to the conclusion that the elements Cy of C_ have the 
form of a radial integral of /(r) multiplied by a product 
Fi{r)Fj{r) Since the (triplet) F-functions arc al- 

most equal over the range where (y**'^)''° contributes sig- 
nificantly, the matrix representation of the spin-angular 
part (fT3|) is the only part in which the dependence on i, j 
survives. We conclude that we can completely leave the 
detailed radial form of (1/^^)''° out of consideration and 
specify only a single multiplicative constant C, its ef- 
fective strength, to be compared to CMies, the strength 
following from the ah initio calculation by Mies et al. [17| . 

An application of this approach to Cs atoms is shown 
in Fig. [5] For atoms lighter than Rb, the magnetic dipole- 
dipole interaction (V^^)'^ adequately describes the de- 
cay of ultracold spin-stretched alkali gases. The second- 
order spin-orbit interaction (y^*')'^°, however, increases 
even more strongly with Z than the first-order spin-orbit 
coupling and is expected to dominate the influence of 
^^ss-j^ for cesium. The figure shows the experimental 
rate coefficient G for the decay of the spin-stretched state 
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(I/, m/ = 4,4)), measured by Soding et al. [23 ■ It also 
shows our calculated rate coefficient as a function of v^t 
for several values of C/CMies- The theoretical range for 
vdt is taken from [l8| . Clearly, the magnetic dipole in- 
teraction on its own (C = 0) is far too weak to account 
for the measured decay rate. The strength constant C of 
j^yss^so j^g^g |-|g factor of about 4 larger than the ab 
initio value CMies to obtain agreement with the experi- 
mental rates [lllli. 



(A) 



III. ADIABATIC ACCUMULATED PHASE 
METHOD 

The theoretical precision needed for the "state of the 
art" BEC and Fermi degeneracy experiments forces us 
to shift ro to atomic distances significantly smaller than 
16ao to neglect the singlet-triplet coupling for r < tq 
according to the above-mentioned condition (1) for the 
applicability of the straightforward accumulated phase 
method. We then run a real risk of violating condition 
(2), however. In this section we present a more sophisti- 
cated variant of the accumulated phase method, already 
introduced briefly in Ref. [llj, that allows us to relax 
condition (1) to some extent, making it possible to find 
a value for 7*0 while achieving the desired accuracy. 

In Fig. [5] we explain the difference between the conven- 
tional accumulated phase method and the new approach, 
distinguishing several intervals along the r axis according 
to the relative magnitudes of y'^^ and V^^"^^. In part A 
we consider three intervals illustrating the conventional 
method. In the left interval V^^ is so weak compared to 
V'"'^, i.e., to the 5 = ^ = 1 splitting of potential 
curves, that the coupling due to V^^^ can be neglected. 
We thus have S = and 1 as a good quantum number. 
The remaining part V^^'^ , together with the two-atom 
Zeeman interaction V"^, can therefore be included effec- 
tively in the Hamiltonian via its eigenvalues, that can 
simply be added to the singlet and triplet potentials, in 
addition to their centrifugal I splitting. We thus have a 
set of singlet and a set of triplet potential curves, each 
with known energy separations independent of r. In the 
right interval of part A the situation with respect to the 
relative magnitude of V^^ and V''^'^^ is opposite and the 
individual atomic hyperfine labels /i, /2, '7i/2 char- 
acterize the spin states. In the middle interval the two 
potential terms are comparable. The separation radius 
ro is chosen as far right as possible in the V^^ <C V°^^^ 
interval. The boundary conditions for the pure singlet 
and triplet radial wave functions at ro along the poten- 
tial curves can therefore be formulated simply in terms of 
E and I dependent pure singlet and triplet phases (j){E, I) 
(note that E is the energy relative to dissociation). 

The new insight leading to our alternative approach 
concerns the role of V^^~ . Let us turn to part B of Fig. [6] 
and consider what happens when we move into the region 
where V^^ ~ ^oxch^ Qj^^ ^-^^^ ^^^.^^ pass through an in- 
terval where the V^^~ coupling is not negligible but still 



FIG. 6: Subdivision of radial ranges to illustrate choices of ro. 
Part A distinguishes three radial ranges. In the left interval S 
is a good quantum number. In the right interval the individual 
atomic hyperfine labels /i, my i, f2,rnf2 characterize the spin 
states. Conventionally, ro is chosen as far right as possible in 
the V^^^ <^ ycxch ijj^g].ya,I. Part B shows the radial intervals as 
they occur in the adiabatic accumulated phase method. The 
intermediate radial interval is subdivided in one in which the 
influence of V^^~ is small and adiabatic and one in which it is 
not. The separation radius ro is chosen as far right as possible 
in the former interval. 



small and adiabatic. In principle, V induces both 
a spin mixing between the S = and 1 states, and a 
perturbation on the radial wave functions. We include 
the spin mixing, but neglect the radial perturbation, so 
that the radial functions are still decoupled singlet and 
triplet waves characterized by pure singlet and triplet ac- 
cumulated phases. Note that the spin mixing is a first 
order perturbation, whereas the energy perturbation on 
the singlet and triplet states is a second order effect. 

As a further illustration of the difference between the 
two methods we discuss the example of ^^Rb 4- ^^Rb 
scattering with initial spin state /2, "^/2 > = 

|1, -1, 1,-1 >. We calculate the matrix £ of Eq. (gS]), 
but instead of taking it to arise from the spin-spin inter- 
actions in the range r < ro, we evaluate how it builds up 
from V^^~ by solving a coupled channels equation based 
on y^"" + V^'^^ -I- V'^ for this radial regime. We note that 
in the present application of Eq. ((23|) the F and G waves 
continue to arise from the part 1/^^°" + -|- ]/Z q£ ^j^g 
Hamiltonian. 

The solid line in Fig. [7] shows the largest C matrix el- 
ement in absolute value for rg values in the range [11.65, 
16.0] oo from the latter calculation. The dashed line is 
the analogous quantity following from the adiabatic accu- 
mulated phase method. Clearly, the latter is in excellent 
agreement with the rigorous result for the small rg val- 
ues. The error gradually grows to 0.25 x 10~^ at ro = 
16.0 qq. This amounts to an error of about 10 % of the 
total effect due to V^^~ , which by itself is of order 0.4 % 
of the analogous V^^+ quantity (j)^ x E^f {^"^ Rb) ~ 0.6. 
Note that the conventional accumulated phase method 
corresponds to Cij = 0. 

We emphasize that the new approach includes effec- 
tively the adiabatic spin mixing in the complete range 
r < rp. Although we impose the boundary condition 
that starts the coupled channel calculation in the range 
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FIG. 7: Comparison of conventional accumulated phase 
method and new alternative approach for *^Rb + *'^Rb 
scattering with initial spin state m/i, /2, m/2 > ~ 
|1,— 1,1,— 1 >. Solid line: largest C matrix element dj 
in absolute value for full coupled channel calculation in ro 
range [11.65, 16.0] oq. Dashed line; analogous result from 
adiabatic accumulated phase method. Conventional method 
corresponds to dj = 



r > To only at ro, by its local character the adiabatic spin 
mixing may be understood to have been included for all 
smaller r values. This is clearly illustrated in Sec. IVIl 
where we discuss an application of the adiabatic accu- 
mulated phase method to ^^Rb and ^^Rb, previously pre- 
sented in Ref . . It turns out (see column C of table 
H] in the following) that the deduced potential parame- 
ters and vds^vdt are highly independent of tq over a 
rather long range. An important aspect is a comparison 
with the straightforward accumulated phase method. In 
particular, we will present convincing evidence, in addi- 
tion to Fig. [71 that the new variant allows us to shift tq 
to larger interatomic distances without significant loss of 
accuracy, thus enabling us to use more reliable potential 
terms in the range of interatomic distances r > tq in the 
form of dispersion and exchange expressions with a small 
number of parameters. 



IV. MASS SCALING: EXPLICIT ISOTOPIC 
DEPENDENCE OF PHASE PARAMETERS 

As long as experimental data are analyzed for bound 
states and cold collisions of a single pair of (un)like atoms, 
it is only the local phase at tq, i.e., the modulo tt part of 
the accumulated phase I) that plays a role in the ra- 
dial boundary condition. In this section we consider the 
combined analysis of several isotopic versions of atom 
pairs and the advantages of mass scaling in that connec- 
tion. We believe that this subject will play an increas- 
ingly important role in cold atom physics, also for colli- 



sions of unlike atoms [29| . When analyzing experimental 
data for two isotopic pairs, making use of the first terms 
of the Taylor expansion (|17p , we would need to introduce 
a set of 2 (5 = 0, 1) times 3 (0°, (p^, and (/)') indepen- 
dent parameters for each of the two-atom systems, to be 
determined by comparing theoretically predicted to ex- 
perimentally determined properties of cold collisions or 
weakly bound states. 

The mass scaling is based on both the Born- 
Oppenheimer and WKB approximations. The former 
approximation enables us to assume equal central poten- 
tials for the isotopic pairs. Clearly, it is essential for this 
approach that Born-Oppenheimer breakdown corrections 
can be neglected. The WKB approximation makes it pos- 
sible to use an explicit expression for the accumulated 
phases as radial integrals containing the reduced mass 
via the wave number fc(r). As we will see, the actual 
value for ro chosen in applications of the adiabatic ac- 
cumulated phase method is at small enough interatomic 
distances along the outer slope of the potentials wells 
for the relative atomic motion to provide for an accurate 
validity of the WKB approximation in the radial range 
r < rg. We start from the WKB integral above, 
written more specifically as 



HE, I) 



k (r) dr 



(24) 



with the inner turning point and the added constant 
7r/4, associated with the quantum mechanical penetra- 
tion into the inner wall of the potential [2^(Ch. VI). We 
thus have the proportionalities 



5 « 



(25) 



and by differentiation of Eq. 



with respect to E and 



1,1 _ 



dE 



d(j) 



1=0 



dl{l + l) 



E=0 



fidr 



dr 1 

OC . 

2fcr2 ^ 



(26) 



(27) 



Clearly, the advantages of a combined analysis of iso- 
topes and the associated mass scaling are a) We extend 
the set of available experimental data without increasing 
the number of fit parameters: we need the phase param- 
eters of only one of the isotope pairs, b) Via the scaling 
of the fit becomes sensitive to the number of nodes 
of the radial wave function left of ro, in addition to the 
modulo TT part of (f)^. With the dispersion + exchange 
parameters deduced in the analysis we then also know 
the number of nodes on the right-hand side and thus the 
numbers of bound singlet and triplet states for all pos- 
sible isotope pairs, not only those analyzed. We will see 
an example of this approach in the case of ^^Rb-f*'^Rb 
and 87Rb-h87j^i3 |vn] 
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Eqs. ((26l) and (jST)) enable us to mass-scale and for 
two isotopic pairs A = A-y.A^ and yl' = A!^ , (A, , A!^ 
standing for atomic mass numbers) according to 



n 



and 



-1 A' j\ 



(28) 



where TZ = \flKAjiKAJ with [i. being a reduced mass. For 
these scaling equations contributions to 0(-E, V) indepen- 
dent of E and I do not play a role. For the mass scaling 
of 0*^, on the other hand, we have 

= + (29) 

with rif,' the number of zero-energy s-wave nodes up to 
the radius of interest (ro), excluding the node at r = 0, 
and the modulo tt part of the total phase 0°. 

Each phase cycle tt corresponds to one additional radial 
node and thus an extra (vibrational) bound state in the 
potential. 

Combining this equation with the mass scaling relation 
(I25l) we find 



UhTT =TZ 



Vmod(TT) 



A' I 

^ 4 



(30) 

so that the scaled ^^q^^j^) values of the two isotopic pairs 
are related according to 



n ■ 



moa{7T) 



(31) 

and its inverse, obtained by interchanging the isotopic 
atom pairs and substituting 1/7?. for TZ. The last term 
gives rise to a number of discrete values for the mass- 
scaled modulo TT phase of isotopic atom pair A, depending 
on nj, for the other pair. The interval between these 
discrete values is (1 — TZ)tt. This discretization can be 
exploited when extracting information from experimental 
data of multiple isotopic pairs and requiring the modulo 
TT phases for the pairs considered to be related according 
to Eq. ((3T|) . Clearly, this allows us to deduce n'^ and, 
by exchanging the roles of the isotope pairs, "^nj^. It 
should be emphasized that the (adiabatic) accumulated 
phase method thus offers a unique possibility to deduce 
numbers of bound states for potentials without knowing 
their short-ran ae p art up to tq. This approach has been 
applied in Ref. in the analysis of a set of experimental 
^^Rb and ^^Rb bound state and cold collision data. In the 
present paper we build on that analysis, which we wish 
to describe and discuss in more detail. We come back 
to this in connection with column A of Table [I] that has 
been taken from [ll[. In the same context we estimate 
the accuracy of the mass scaling for these isotopes. 



V. ACCURACY OF MASS SCALING 

A crucial issue for the possibility to combine the anal- 
ysis of different isotope pairs is its expected accuracy. 
In that connection two types of corrections need discus- 
sion, corresponding to the adopted Born-Oppenheimer 
and WKB approximations. 



Accuracy of mass scaling: adiabatic correction to BO 

The main correction to the Born-Oppenheimer approx- 
imation is the adiabatic or diagonal correction Vad to the 
interatomic potential (30j , given by 

1 

Vad{r) =< i^ei{x]r)\ - —l\r\i}^i[x]r) >cx -, (32) 
2/x /X 

with ipei the electronic wave function (x = electronic co- 
ordinates), depending parametrically on the nuclear co- 
ordinates. This leads to an adiabatic correction to the 
accumulated phase 



4'adiE, I) = --2 



dr 
Mr) 



Vadir). 



(33) 



To show its classical meaning we write it as a time in- 
tegral over the collision in the classically allowed range 
within rg {dt = dr/v{r[t])): 



^adiEJ) 



Vadir[t])dt 



1 



ol < Vad >cl, 



(34) 

proportional to 1/ The last member of this equation 
indicates the proportionality to the collision time TcoI and 
to a classical expectation value in this range. In the fol- 
lowing we estimate the isotopic spread ISVad and thus 
the associated spread in accumulated phase parameters 
on the basis of experiment, on the basis of theory, and 
using a combination of both. 
1. Experimental evidence 

In 2000 a paper by Seto et al. [3l| described a measure- 
ment of high-resolution A — s- A" emission data for a mix- 
ture of the isotopic pairs *^5Rb2, ®'^Rb2, and ^^Rb^^Rb, 
covering in total 12148 transition frequencies. The data 
allowed a ground-breaking analysis of vibrational level 
spacings of the A^E+ electronic state up to w = 113 (r 
up to 25ao). Although the data set, with uncertainties 
±0.001 cm^^, involved the above three isotopic pairs, the 
analysis turned out to lead to a common singlet potential 
without any sign of a Born-Oppenheimer breakdown. A 
similar analysis for the triplet case does not exist. 

This result enables us to deduce an upper limit for the 
correction to a mass scaled singlet phase due to Born- 
Oppenheimer breakdown. To that end we consider the 
isotopic difference A(f>ad{E, I) of the adiabatic phase cor- 
rection and note that the above ±0.001 cm~^ uncertain- 
ties correspond to quantum mechanical expectation val- 
ues of the isotopic difference AVadir) over a large set of 
rovibrational states v, I with probability densities cover- 
ing together at least the whole range [rt, ro]. This justifies 
the conclusion that the isotopic difference AVadir) is less 
than 0.001 cm^^ in absolute value. For energies E close 
to and using Eq. we thus find a correction due to 
the implicit isotopic dependence: 



\A(j)°s\ < 0.001cm' 



if 



0.33 X 10" 



(35) 
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Here and in the following these estimates apply to the 
isotopic pairs ®^'®^Rb2 - ^^'^^Rb and half these values to 
the pairs ss-ss^b^ _ s^^si^h and ^s.s^j^ba - ^^-sr^b. We 
have used the value of d(t>s/dE = from the analysis 
in Ref . [ll[ . In the final result we have split off a factor 
TT representing the basic periodicity associated with the 
phases (p. We expect a similar order of magnitude for the 
implicit isotopic correction in the triplet case. 
2. Theoretical evidence 

An order of magnitude estimate for both the singlet 
and triplet case can be based on the long-range exp ression 
for Vad proposed by Dalgarno and McCarroll [33 |: 



Vad = --r- 



Vbo + r 



dVBo{r) 
dr 



(36) 



with nie the electron mass and Vbo = Vs/t the Born- 
Oppenheimer potential for the atom pair. Assuming that 
P6p can be used for an order of magnitude estimate in 
the range [rt , ro] [s^l , we thus obtain 



A4>ad{E,l) 



rrie A/i 



1 



k{r) 



V{r) 



dV{r) 



dr. 



With the singlet potential of Ref. 
triplet potential from Ref. 



dr 

\ (37) 

3ll | and the ab-initio 
3j| for r < ro, both shifted 



'vertically' and smoothly joined to dispersion ± exchange 
forms following from the parameters in the later Table U 
for r > To, we find 



A, 



<?!)^ = +0.037 X lO^^TT, 
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Av(l>T = -0.19 X 10 -TT. 

(38) 

We note that the smallness of the estimated singlet phase 
correction is due to the large negative contributions to 
the radial integral over the Dalgarno-McCarroU expres- 
sion (j36p at small r values, which compensate the positive 
contributions at longer range to a considerable extent. 
3. Combined evidence 

To improve the above estimates on the basis of experi- 
ment and theory together, we note for the singlet case 
that the Dalgarno-McCarroU expression ([36]) is larger 
than the maximum adiabatic correction 0.001 cm~^ in 
absolute value allowed by experiment [sH in a range of 
atomic distances starting from the inner classical turning 
point rt = 5.9 oq until 7.7 ag. We therefore use the exper- 
imental limit in the radial integral (j37p until a distance 
of 7.7 oo so that it fits continuously to the theoretical 
prediction in the further interval up to the final radius 
To = 16 ciQ. For the triplet situation rt is much larger 
(about 9.5 ao). In that case the Dalgarno-McCarroU ex- 
pression is smaller in absolute value than 0.001 cm~^ over 
the whole interval [rt, ro]. Substituting that in the radial 
integral, we find our triplet result. In total we find 

|A(/)^| = 0.61 X 10~^7r, |A0?.| 0.19 X IQ-^TT. (39) 



Accuracy of mass scaling: corrections to WKB 

The order of magnitude of this correction is easily es- 
timated by comparing the mass scaled ^^Rb phase pa- 



TABLE I: Interaction parameters (au) derived from combined 
^'Rb and *^Rb experiments (column A) including error bars, 
mainly due to 10% uncertainty in Cio; column B: fractional 
changes due to phase corrections; column C: percentages of 
variation of same quantities over range [10.85,16] ao of ro val- 
ues according to adiabatic accumulated phase method; col- 
umn D: same for conventional method. 



Quantity 

Ce/lO^ 

Cs/lO^ 

Cio/10" 

J. 10^ 
aT(**^Rb) 
os(*'Rb) 

OT(*^Rb) 

as C'Rb) 
«DT(mod l),nbT(*''Rb) 
UDs(mod l),n6s(*''Rb) 
UDT(mod l),n(,T(*^Rb) 
i;Dg(mod l),7i65(^'^Rb) 



A 



B(%) C(%) D(%) 



4.703(9) 
5.79(49) 
7.665 (Ref. [35]) 
0.45(6) 
+98.98(4) 
+90.4(2) 
-388(3) 
+2795lg° 
0.4215(3), 41 
0.455(1), 125 
0.9471(2), 40 
0.009(1), 124 
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rametcrs to those obtained by numerical integration of 
the singlet and triplet radial Schrodingcr equations up 
to r = 16ao for the above-mentioned singlet and triplet 
potentials with the reduced masses involved. The devia- 
tions of the mass scaled phases are 



\A(j)°s\ = \A(j)'^\ = 2 X lO-^TT. 



(40) 



Comparison of phase corrections to error bars from 
previous analysis 

To illustrate the smallness of the above estimated 
phase corrections, we compare them with the error bars 
obtained in our previous brief descri ptio n of the adiabatic 
accumulated phase method in Ref. [ll| . In that Letter a 
combined analysis of *^Rb and ^^Rb experimental data 
led to values for interaction and scattering properties of 
Rb atoms with an unprecedented accuracy. In column A 
of Table [J we recapitulate the dispersion coefficients Cg , 
Cg, the strength parameter J of the exchange interaction, 
and the set of pure singlet and triplet scattering lengths 
+ associated fractional vibrational quantum numbers at 
dissociation vjj, together with their error bars. Column 
B gives for comparison the maximum fractional changes 
(in %) of the same quantities that result from the com- 
bination of the two types of phase corrections above. We 
conclude that the latter are small compared to the error 
bars resulting from the analysis in Ref. [ll| and indi- 
cated in column A. The latter are mainly due to the 10 
% error assumed for the theoretical Cio value taken from 
Ref. [11]. The largest of the fractional phase corrections 
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is that for J. We note that that is not unexpected tak- 
ing into account that this concerns the coefficient of a 
radiaUy exponential term, which is extremely sensitive 
to the damping coefficient in the exponential. This also 
explains the relatively large error bar in column A. 

The beautiful agreement with experiment, achieved 
in the analysis of Ref. [ll| . is a convincing further in- 
dication that the mass scaling procedure is an excel- 
lent approximation. For instance, the values of Cg and 
Cg agree with values Cg = 4.691(23) x 10^ [sl and 
Cs = 5.77(8) X 10^ [131, calculated by Derevianko and 
co-workers, while J agrees with the most recent calcu- 
lated value J = 0.384 x 10^ published by Hadinger and 
Hadinger [3]. 

We can also conclude that there is considerable room 
for an extension of the mass scaling procedure to appli- 
cations of the adiabatic accumulated phase method to 
isotopic pairs of lighter elements than the Rb isotopes 
studied here, despite the expected larger phase correc- 
tions due to Born-Oppcnhcimcr and WKB breakdown. 



VI. COMPARISON OF ADIABATIC TO 
CONVENTIONAL ACCUMULATED PHASE 
METHOD AND DEPENDENCE ON ro 

To illustrate the advantages of our adiabatic accumu- 
lated phase method, we compare a calculation including 
the adiabatic spin mixing at tq to one without, i.e., the 
conventional approach. In both cases we consider the 
optimization of the potential parameters given a set of 
^^Rb and ^^Rb experimental data according to the anal- 
ysis in Ref. [ll|. It turns out that the optimized values 
of the quantities in Table |T] are highly independent of the 
choice of tq. To demonstrate that we have added the 
percentages of variation over the tq interval [10.85,16.0] 
in column C of the table. In this case too the exchange 
strength parameter J is an exception, with a variation 
of 1%. This can be explained as indicated above in con- 
nection with columns A and B of Table [J In column 
D we have added for comparison the significantly larger 
percentages of variation of the same quantities according 
to the conventional accumulated phase method. 

In Fig. [8] we select by far the most sensitive param- 
eter, J, to show its To dependence. Each of the curves 
shows the behavior of the fractional deviation AJ/Jth 
of the optimal J value from the value Jth published by 
Hadinger et al. [13] ■ The + signs connected by the dashed 
curve show the result of a calculation along conventional 
lines. Each point indicated on the curve represents the 
outcome of a separate optimization. Switching on the 
spin mixing adiabatically at tq gives rise to the solid line. 
Clearly, the oscillation is strongly reduced. The remain- 
ing oscillation is mainly due to the WKB correction and 
the non-adiabaticity of switching on the coupling due to 
V^^^ . Even shifting ro to 16ao keeps the oscillation am- 
plitude in J to below the 1% level. 

Fig. [5] suggests that one might just as well select a 
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FIG. 8: Fractional correction to the literature value Jth for the 
strength parameter of the exchange energy versus ro [a. The 
dashed line connects points calculated with the traditional 
accumulated phase method, the solid curve to those using the 
adiabatic accumulated phase method. 

smaller value for ro near 12 ao to avoid the V'^^" cou- 
pling issue altogether. If we would have done that from 
the beginning, however, we would have missed a key mes- 
sage from our study: the fact that the final results are 
highly independent of the central potentials within an 
interatomic distance of 16ao. This applies in particular 
to the exchange potential V'^^'^^ for which the Smirnov- 
Chibisov radial dependence ([6]) is an asymptotic expres- 
sion. The same applies to the asymptotic expression ([S]) 
for the dispersion potential. 



VII. DETERMINING NUMBERS OF SINGLET 
AND TRIPLET BOUND STATES FOR **^RB + 
**^RB AND *^RB + *^RB SYSTEMS 

Here we come back to the relation between the mass- 
scaled modulo TT accumulated phases for different iso- 
topic versions of a general atom-atom system discussed 
in Sec. IIVI in particular Eq. |3T|) . This relation and its 
inverse contain the (unknown) numbers of nodes n'^ of the 
zero-energy radial wave function contained in the poten- 
tial from the inner turning point up to rg for the two 
interrelated atom pairs A and A! . As pointed out above, 
this enables us to deduce the total numbers of bound sin- 
glet and triplet states from available experimental data. 

To illustrate our method via the example of the ^'^Rb 
+ *^Rb and ^'^Rb -I- ®''Rb systems, it is helpful to reca- 
pitulate some aspects of the analysis that led to column 
A of Table [J in Ref. [IH , reproduced in the above Table 
m The experimental material analyzed consisted of data 
on cold collisions and on bound states exceptionally close 
to the continuum, partly for ^^Rb and partly for ^^Rb. 
The six parameters varied in a analysis were ^^^^j^. 
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, ^^(/)^, Ce, Cs, and J, with Cio held fixed at the 
theoretical value from Marinescu et al. [s^ . 

Instead, in part A of Fig. [9] we present a contour plot of 
the reduced value for the subset of ^^Rb data as a func- 
tion of Ce and ^^DDT(niod 1) [equivalent to '^^^^.(mod 
tt)], keeping the remaining parameters , ^''^^p, Cs, 
and J fixed at their values according to column A. The 
dispersion coefficient Cg is expected to be within the in- 
terval indicated by the horizontal double arrow [sBl ■ The 
absolute minimum is the white point indicated within 
the < 10^ area. Part B is based on the ^^Rb data only, 
again keeping , *^c/>?p, Cs, and J fixed and calcu- 
lating the three ^^Rb phase parameters by mass scaling. 
The dashed line shows the bottom of the 'deepest trench' 
in the surface of part A, with the overall minimum in- 
dicated by the square. The solid lines indicate equivalent 
trenches in the surface for the ^^Rb data. Note that 
the interval between the discrete *^Rb phase values in 
Eq. dSI]) is approximately (1 - 7^)7r w 0.0127r w 0.036, 
corresponding to intervals (1 — 7?.) w 0.012 for ^''vot 
between the solid lines. The position of the minimum 
X^ value along each trench is indicated by a solid circle. 
Each of the solid lines corresponds to a specific number 
^''n'f, of nodes assumed to be contained in the ^''Rb triplet 
potential up to ro = 16ao: ^"^nj, 19,. ...,25. Clearly, 22 
is the preferred value for ^^nj,. It turns out that the lower 
mass in the ^'^Rb case does not decrease this value: 
~ 22. With the optimalized dispersion and exchange 
parameters of Table U we can continue the solution of 
the zero energy triplet radial Schrodinger equations for 
the two atom pairs to find the total numbers of bound 
states. For the case of ^^Rb we thus find the total num- 
ber of triplet bound states nbT(^^Rb2) = 41. Similarly 
for s^Rb: = 40. With the singlet potential of 

Seto et al., suitably shifted vertically and also completed 
with a long range part the result is ~ 125, 

n;,s(''Rb2) = 124. 



VIII. SUMMARY AND OUTLOOK 

We have presented a theoretical method that enables 
one to describe and predict the interaction and scattering 
properties of (ultracold) atoms. It allow us, for instance, 
to predict the ®^Rb spinor condensate to be ferromag- 
netic [m, a prediction for which the relevant scattering 
lengths have to be calculated with a precision better than 
1%. It is also unparalleled in comprehensiveness: it al- 
lows the prediction of a large and varied set of experi- 
mental quantities for all pairs of like and unlike atoms. 

Its original version, the accumulated phase method, 
was designed to predict essential properties like scatter- 
ing lengths and Feshbach resonances, enabling the real- 
ization of Bose-Einstein condensates and Fermi degen- 
erate gases of alkali atoms, for which the short-range 
interaction was insufficiently known to calculate these 
properties directly. The method consisted of replacing 
the short-range interaction with a boundary condition 
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FIG. 9: (A) Contour plot of the reduced as a function 
of Ce and *''^?;_DT(niod 1) for *^Rb data only. The dispersion 
coefficient Ce is expected to be within the interval indicated 
by the horizontal double arrow [s^. (B) The dashed line in- 
dicates the bottom of de 'deepest trench' in the surface of 
part (A), with the overall minimum indicated by the square. 
The solid lines indicate equivalent trenches in the x^ surface 
for *^Rb data only, making use of mass scaling to deduce 
*^(^y(mod tt) from *'^i'_DT(mod 1). The position of the mini- 
mum x^ value along each trench is indicated by a solid circle. 
Each of the solid lines corresponds to a specific number *^nj, 
of nodes assumed to be contained in the ^'^Rb triplet poten- 
tial up to ro: *^nj, — 19,.... ,25. Both graphs A and B have 
been generated by varying *^DDT(mod 1) and Ce, keeping 
the remaining fit parameters *^(jf>^, '^'^ <f)\ J, and Cg of the x^ 
analysis in Ref. [ll[ fixed at their optimal values. 



on the two-atom wave function at an interatomic dis- 
tance r = ro, deducing the boundary condition from 
available experimental data, and predicting all other rel- 
evant data. The new, adiabatic, version of the method, 
described in the present paper, has been presented briefiy 
in a previous Letter [ll[. Whereas the original method ne- 
glected the hyperfine coupling between singlet and triplet 
states for r < ro and included this coupling together with 
asymptotic dispersion -I- exchange expressions for r > tq , 
the new approach takes the adiabatic singlet-triplet mix- 
ing by Vhf into account at the separation radius tq and 
therefore effectively also at smaller r, neglecting the (sec- 
ond order) changes of the radial waves. This makes it 
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possible to shift rg to larger interatomic distances, thus 
allowing for more reliable asymptotic potential terms in 
the range r > rg. 

We have described how the second-order spin-orbit 
interaction can be included, as well as a mass scaling 
approach to relate the accumulated phases for differ- 
ent isotopic versions of atom pairs. The accuracy of 
the mass scaling has been discussed, taking into account 
both Born-Oppenheimer and WKB breakdown. Esti- 
mates have been given for the Rb isotopes, pointing to 
a high accuracy. Using the Rb isotopes for illustration, 
the adiabatic and conventional accumulated phase meth- 
ods were compared, and the tq dependence of their op- 
timized interaction parameters was studied. Finally, we 
have explained how the total numbers of bound singlet 
and triplet two-atom states follow from a combined anal- 
ysis of different isotopic versions of atom pairs, without 
knowing the short range interatomic interaction. 

We believe that the new adiabatic accumulated phase 



method has great potential to support further studies 
of cold atom systems, especially in the rapidly grow- 
ing field of pairs of unlike atoms, to which the method 
can readily be extended [2^. The set of phase param- 
eters that it makes use of can be systematically ex- 
tended when larger energy or angular momentum ranges 
come into play experimentally, contrary to other choices 
used for the adjustment of the short-range part of model 
potentialsfUlij,!!!. 
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